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ABSTRACT 

In  this  paper  we  study  several  problems  concerning  the  visi- 
bility of  a  polyhedral  terrain  cr  from  a  point  (or  several  points) 
lying  above  it.  Our  results  are:  (1)  For  a  fixed  viewing  point  a, 
one  can  preprocess  cr  in  time  0{noi{n)  log  n),  to  produce  a  data 
structure  of  size  0(na(n)  log  n),  which  supports  fast  ray  shoot- 
ing queries,  where  each  such  query  asks  for  the  point  on  ct  that 
is  visible  from  a  in  a  specified  direction.  Here  n  is  the  number 
of  faces  of  u  and  a(n)  is  the  extremely  slowly  growing  func- 
tional inverse  of  Ackermann's  function.  (2)  If  the  viewing  point 
a  can  vary  along  a  fixed  vertical  line  L,  then  the  entire  visibility 
structure  of  cr  from  L  is  of  combinatorial  complexity  0{n\<i{n)), 
where  \n(n)  is  the  maximal  length  of  an  (n,4)  Davenport- 
Schinzel  sequence,  and  is  nearly  linear  in  n,  and  where  the  visi- 
bility structure  in  question  is  the  decomposition  of  Lx5^  into 
maximal  connected  regions,  such  that  for  each  such  region  R,  all 
points  (a,u)  €  R  are  such  that  the  ray  from  a  €  L  in  direction 
M  6  5^  first  intersects  a  at  a  point  on  the  same  face  of  a.  Furth- 
ermore, we  present  an  C>(nX4(n)  log  n)-time  algorithm  that 
preprocesses  L  and  cr  into  a  data-structure  of  size  0(n\4{n)) 
which  supports  O(log^n)  time  ray  shooting  queries.  (3)  Con- 
cerning the  results  in  (2)  we  show  that  (i)  if  L  is  not  vertical, 
then  the  resulting  visibility  structure  can  be  of  size  Cl{n^);  (ii) 
there  exist  a  vertical  line  L  and  a  polyhedral  terrain  cr  with  n 


Work  on  this  paper  has  been  supported  by  Office  of  Naval  Research  Grant  N00014-82-K-0381,  Nation- 
al Science  Foundation  Grant  No.  NSF-DCR-83-20085,  and  by  grants  from  the  Digital  Equipment  Cor- 
poration, and  the  IBM  Corporation.  Work  by  the  first  author  has  also  been  supported  by  NSF  Grant 
DCR-84-01633  and  by  an  IBM  faculty  development  award. 


-2- 


faces,  for  which  the  resulting  visibility  structure  is  of  size 
n(n^a(n)).  (4)  Finally  we  consider  the  problem  of  placing  on 
the  surface  a  one  or  several  viewing  points  which  collectively 
cover  the  entire  surface  (i.e.  each  point  on  a  is  visible  from  at 
least  one  of  these  viewing  "stations").  We  show  (i)  in  the  case 
of  a  single  viewing  station,  one  can  determine  in  time 
0{n  log  n)  whether  such  a  station  exists,  and  if  so  produce  such 
a  point;  (ii)  the  problem  of  finding  the  smallest  number  of 
points  on  a  that  can  collectively  see  the  entire  surface  cr  is  A^P- 
hard. 


1.  Introduction 

In  this  paper  we  consider  the  variety  of  3-D  visibility  problems  men- 
tioned in  the  abstract,  and  develop  efficient  techniques  for  their  solutions,  as 
outlined  above.  These  problems  can  be  viewed  as  simple  3-dimensional  gen- 
eralizations of  planar  visibility  problems.  For  example,  Chazelle  and  Guibas 
[CG2]  (cf .  also  [GHLST])  consider  several  ray  shooting  problems  for  points 
lying  inside  a  simple  polygon  P.  The  results  in  [CG2]  and  in  [GHLST]  show 
that,  after  triangulation,  P  can  be  preprocessed  in  linear  time  to  support 
logarithmic-time  ray  shooting  queries  in  which  the  point  from  which  the  ray 
emerges,  as  well  as  the  direction  of  the  ray,  can  vary  from  one  query  to 
another. 

A  direct  generalization  of  this  result  would  aim  to  preprocess  a  simple 
polyhedron  cr  so  as  to  support  fast  ray  shooting  queries  from  arbitrary  inte- 
rior points  in  arbitrary  directions  (these  are  queries  in  which  the  first  inter- 
section of  such  a  ray  with  the  boundary  of  ct  is  sought).  This  problem  how- 
ever seems  to  be  much  harder  than  the  two-dimensional  case  if  efficient 
preprocessing  time  and  space  are  desired.  To  obtain  a  preprocessing  pro- 
cedure  whose    complexity    resembles    that    of    the    solution    in    the    two- 


dimensional  case,  the  problem  has  to  be  restricted  considerably.  The  first 
result  of  this  paper,  described  in  Section  2,  is  an  efficient  preprocessing  tech- 
nique for  the  case  in  which  cr  is  a  polyhedral  terrain  (i.e.  a  polyhedral  surface 
having  exactly  one  intersection  with  each  vertical  line),  and  the  origin  of  the 
query  rays  is  fixed  at  some  point  a  lying  above  ct.  In  this  case  we  obtain  an 
algorithm  requiring  preprocessing  time  and  space  0{n  a(n)log  n),  which  sup- 
ports C>(log  n)-time  ray  shooting  queries;  here  a{n)  is  the  extremely  slowly 
growing  inverse  Ackermann's  function.  This  is  a  substantial  improvement 
over  a  naive  approach  of  calculating  a  decomposition  of  the  space  of  all  rays 
emerging  from  a  (a  space  conveniently  represented  by  the  unit  sphere  S^) 
into  connected  regions,  each  consisting  of  rays  all  hitting  first  the  same  face 
of  a.  Indeed,  it  is  well  known  that  such  a  decomposition  can  consist  in  the 
worst  case  of  n(n^)  regions  (see  Fig.  1.1  for  a  simple  example  of  this  kind). 

This  result  has  potential  applications  to  computer  graphics  and  is  closely 
related  to  the  hidden  surface  removal  problem  (see  [MK],  [De]).  It  can  be 
used  to  provide  compact  representation  of  an  image  of  a  polyhedral  terrain  cr 
as  viewed  from  some  point,  from  which  the  actual  pixel-by-pixel  image  can 
be  quickly  generated,  or  the  portion  of  cr  appearing  at  certain  "query  pixels" 
can  be  determined  efficiently  without  having  to  generate  the  entire  image. 

When  the  origin  a  of  the  shooting  rays  is  allowed  to  vary,  the  problem 
becomes  harder  (even  if  cr  is  still  assumed  to  be  a  polyhedral  terrain,  rather 
than  an  arbitrary  polyhedron),  and  we  have  not  been  able  to  obtain  a  compar- 
ably efficient  preprocessing  procedure  for  these  extended  problems.  In  Sec- 
tion 3  we  consider  such  an  extension,  in  which  the  point  a  can  vary  along 


some  fixed  vertical  line  L.  Let  us  denote  by  the  visibility  structure  of  a  from 
L  the  decomposition  of  LxS^  into  maximal  connected  regions;  for  each  such 
region  R,  for  all  points  {a,u)  6  /?,  the  ray  from  a  ^  L  in  direction  u  i  S^ 
first  intersects  cr  at  a  point  on  the  same  face  of  ct.  As  stated  in  the  abstract 
we  show  that  the  combinatorial  complexity  of  this  visibiUty  structure  is  at 
most  0(n\^{n)),  where  X4(n)  is  the  maximal  length  of  an  (n,4)  Davenport- 
Schinzel  sequence,  and  is  thus  nearly  linear  in  n.  See  [HS],  [Shi],  [Sh2]  for 
details  about  Davenport-Schinzel  sequences  and  their  analysis. 

We  also  present  an  0(n\4{n)  log  n)-time  algorithm  that  preprocesses  L 
and  a  into  a  data-structure  of  size  0(n\^(n))  which  supports  O(log^n)  time 
ray  shooting  queries. 

To  supplement  the  results  just  mentioned  we  show  that 

(i)  If  L  is  not  vertical,  then  the  resulting  visibility  structure  can  be  of  size 
n(n^).  .,-..-.. 

(ii)  There  exist  a  vertical  line  L  and  a  polyhedral  terrain  cr  with  n  faces,  for 
which  the  visibility  structure  of  Lx 5^  is  of  size  Cl(n^a{n)). 

By  letting  the  viewing  point  a  vary  along  more  general  domains  L,  and 
by  replaving  cr  by  an  arbitrary  polyhedral  scene,  we  obtain  a  collection  of 
progressively  more  complex  versions  of  the  3-D  visibility  problem.  Some  of 
these  generalized  problems,  in  which  the  viewing  point  lies  at  infinity,  are 
discussed  in  [PD]. 

Finally  we  consider  in  Section  4  the  problem  of  placing  one  or  several 
viewing  points  on  the  surface  u,  which  collectively  cover  the  entire  surface 


(i.e.  each  point  on  o-  is  visible  from  at  least  one  of  these  viewing  "stations"). 
Such  visibility  problems  apse  in  many  applications,  and  can  be  thought  of  as 
3-D  generalizations  of  the  well  known  problems  of  stationing  guards  in  art 
galleries  [Ch];  see  pFPAN]  for  a  survey  of  applications  of  this  kind. 

We  show  that 
(i)  In  the  case  of  a  single  viewing  station,  one  can  determine  in  time 
0{n  log  n)  whether  such  a  station  exists,  and  if  so  produce  such  a  point. 

(ii)  The  problem  of  finding  the  smallest  number  of  points  that  can  collec- 
tively see  the  entire  surface  cr  is  NP-haid. 

2.  Visibility  of  a  polyhedral  surface  from  a  fixed  point 

Let  cr  be  a  2-D  polyhedral  surface  in  E^,  that  is  cr  is  the  graph  of  a 
polyhedral  function  z=F(x,y),  defined  over  the  entire  x—y  plane,  and  sup- 
pose cr  has  n  edges  (and  thus  0{n)  vertices  and  faces).  Let  a  =  {xo,yo,zo)  be 
a  fixed  point  lying  above  cr  (i.e.  ZQ>F(xo,yo)).  Our  goal  is  to  preprocess  ct 
into  a  data-structure  which  supports  fast  responses  to  ray-shooting  queries, 
which  ask  for  the  determination  of  the  first  point  on  a  hit  by  a  query  ray 
emerging  from  a  in  a  specified  direction. 

Note  that  each  such  ray  r  can  be  represented  by  a  point  F  =  (0,(t))  on  the 
unit  sphere  5^,  where  4)  is  the  angle  between  r  and  the  positive  z-direction 
(the  "azimuth"  of  r),  and  6  is  the  planar  orientation  of  the  horizontal  projec- 
tion of  r.  Our  problem  thus  calls  for  calculating  a  partitioning  n  of  5^  such 
that  for  each  region  R  of  IT,  the  first  intersection  with  a  by  all  rays  in  R 
emerging  from  a  is  at  the  same  face  of  cr.  Once  IT  is  available,  we  can  further 
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preprocess  it  using  any  one  of  the  known  algorithms  for  point  location  (e.g. 
that  of  [DSST])  so  as  to  support  efficient  (i.e.  0(log  n)  time)  point  location 
queries  in  11,  which  would  then  correspond  to  efficient  calculation  of  the  face 
of  cr  first  hit  by  any  query  ray.  However,  as  noted  in  the  introduction  (cf. 
also  [MK]),  the  partitioning  11  might  consist  of  Cl{n^)  regions  (it  always  con- 
sists of  at  most  0{n^)  regions).  Hence,  instead  of  calculating  11  explicitly, 
we  want  to  develop  an  alternative  more  space-efficient  data-structure  that  will 
still  support  ray  shooting  queries  from  a  in  0(log  n)  time. 

Our  solution  proceeds  as  follows.  For  any  geometric  object  u  in  3-space, 
denote  by  «*  its  projection  onto  \h&  x—y  plane.  The  projection  a*  of  ct  can 
be  regarded  as  a  planar  map  whose  vertices,  edges  and  faces  are  projections 
of  corresponding  vertices,  edges  and  faces  of  cr.  We  split  tr*  into  two  disjoint 
parts  o-*,  cr2  such  that  each  point  in  cri  lies  nearer  to  the  projection  a*  of  a 
than  any  point  in  a^-  More  precisely,  define  a  partial  order  on  the  edges  of 
CT*  so  that  for  any  pair  u,  v  of  edges  of  a* ,  u<v  if  there  exists  a  (horizontal) 
ray  emerging  from  a*  which  intersects  both  u  and  v  such  that  its  intersection 
with  u  is  nearer  to  a*  than  its  intersection  with  v  (see  Fig.  2.1(a)).  To  make 
this  into  a  partial  order  we  need  to  ensure  that  there  is  some  orientation 
(8  =  0,  say)  that  no  edge  crosses;  this  is  easily  adiieved  by  replacing  each 
edge  which  does  cross  6  =  0  with  its  two  subsegments  clockwise  and  counter- 
clockwise from  6  =  0.  (If  we  do  not  have  this  property  there  could  be  a  cir- 
cular sequence  of  edges  e i ,  e 2  ,  .  .  •  ,  ^ t ,  such  that  e^  <  c/  + 1 ,  \<i<k,  and 
e^-^ei  (see  Fig.  2.1(b)).  However,  with  this  extra  property  it  is  easy  to 
check  that  "<"  is  indeed  a  partial  order.)    Moreover,  it  is  easy  to  calculate 
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this  partial  order  by  a  standard  ray-sweeping  procedure  in  time  0{n  log  n), 
and  to  complete  this  order  into  a  total  linear  order  in  additional  0{n)  time, 
using  a  well  known  topological  sorting  procedure.  Let  this  linear  order  be 
ei<c2<  •  ■  •  <fn-  Split  the  edges  of  cr'  into  two  subsets 
(tJ  =  {ei,  .  .  .  ,e„,2}  and  a^  =  {e„t2+i,  •  •  •  >e„}.  Next  calculate  the  upper 
rim  of  the  edges  of  a  whose  projections  lie  in  ct'  (call  this  set  of  edges  ui). 
This  is  a  function  4)=/i(9)  defined  over  the  unit  circle  5\  as  follows.  For 
each  6  6  5  \  let  re  be  the  lowest  ray  that  emerges  from  a,  projects  onto  the 
horizontal  ray  at  direction  6  from  a' ,  and  does  not  contain  any  point  lying 
directly  below  any  of  the  edges  in  cti.  h(B)  is  then  defined  to  be  the  azimuth 
<{)  of  re,  i.e.  the  angle  between  r^  and  the  positive  z-direction.  Qearly  /i(0) 
consists  of  a  finite  number  of  smooth  connected  portions,  such  that  for  each 
such  portion  -y,  all  the  rays  r^  in  -y  pass  through  the  same  edge  of  cti.  The 
endpoints  of  these  smooth  portions  are  at  orientations  9  at  which  either  r^ 
passes  through  an  endpoint  of  some  edge  of  cti,  or  re  passes  simultaneously 
through  two  segments  of  cti.  We  will  also  use  the  term  "upper  rim"  to  refer 
to  the  (not  necessarily  connected)  locus  of  points  on  the  edges  of  cti  touched 
by  the  rays  r^,  for  6  €  5^  See  Fig.  2.2  for  an  illustration  of  the  upper  rim 
of  a  collection  of  segments. 

Having  calculated  this  upper  rim,  we  then  proceed  to  apply  recursively 
the  procedure  just  described  to  the  two  subsets  a  J  and  a  I  of  the  set  of  edges 
of  a*. 

The  output  of  this  preprocessing  is  thus  a  collection  of  upper  rims  of 
subsets  of  edges  of  cr,  arranged  in  a  balanced  binary  tree  structure  T,  so  that 
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the  root  of  T  stores  the  upper  rim  of  the  half  of  the  edges  of  ct  that  lie  nearer 
to  a  (according  to  the  order  defined  above),  and  the  left  and  right  subtrees  of 
the  root  contain  the  collection  of  rims  calculated  recursively  for  the  near  half 
and  for  the  far  half  of  the  edges  of  ct  respectively. 

Having  this  structure  available,  actual  ray  shooting  queries  can  be  pro- 
cessed as  follows.  Let  r  be  a  query  ray  emerging  from  a,  and  let  (6,4))  be  its 
spherical  coordinates.  We  perform  a  binary  search  through  T  by  first  com- 
paring r  against  the  topmost  upper  rim  /i(e),  stored  at  the  root  p  of  T;  using 
an  auxiliary  binary  search  over  6 ,  we  locate  the  smooth  portion  of  h  contain- 
ing the  point  (e,/i(e)),  and  then  compare  (in  constant  time)  ^  against  /i(6);  if 
<{)>/i(e)  then  r  lies  above  all  the  edges  of  cr  in  its  nearer  portion  cti,  and  we 
continue  the  search  through  T  at  the  right  child  of  p;  otherwise  r  must  hit 
some  face  of  ct  lying  in  its  nearer  half,  and  we  continue  the  search  through  T 
at  the  left  child  of  p.  When  this  search  is  completed,  we  will  have  found  two 
edges  ei<e2  which  bound  the  same  face /of  cr  and  such  that  the  horizontal 
projection  r*  of  r  intersects  both  e\ ,  ei  and  r  lies  above  ei  (actually  above  all 
edges  of  a  preceding  cj  whose  projections  intersect  r*)  and  below  ei-  Hence 
the  first  intersection  of  r  with  a  lies  in  /,  and  can  now  be  calculated  in  con- 
stant time.  The  whole  search  takes  O(log^n)  time. 

We  can  reduce  the  search  time  to  0(log  n),  with  only  a  constant  increase 
in  the  space  requirement,  as  follows.  We  use  the  fractional  cascading  tech- 
nique as  described  in  [Co],  [CG],  [DSST],  and  elsewhere.  The  basic  idea  is  to 
reduce  to  a  constant  the  search  time  for  going  from  one  node  in  T  to  the 
appropriate  child.   This  is  achieved  by  subdividing  smooth  portions  of  upper 


rims  into  smaller  angular  intervals  so  that  each  interval  of  a  rim  stored  at  an 
internal  node  overlaps  with  at  most  two  intervals  of  the  rims  stored  in  each 
child.  It  then  follows  that  the  total  search  time  can  be  reduced  to  0(log  n). 
Given  the  original  tree  T,  the  modified  tree  T'  can  be  obtained  in  time  linear 
in  the  size  of  T  by  a  single  sweep  from  the  leaves  to  the  root  of  T,  performed 
level  by  level.  This  construction  is  described  in  [Co]. 

The  complexity  of  the  approach  proposed  above  crucially  depends  on  the 
size  of  the  upper  rims  that  it  calculates.  Fortimately,  even  though  the  way  in 
which  the  entire  ct  is  visible  from  a  can  be  of  Cl{n^)  combinatorial  complex- 
ity, upper  rims  have  almost  linear  size,  as  shown  in  the  following 

Proposition  1:  Let  T={ei,  ...  .e^}  be  a  collection  of  n  (non-intersecting) 
segments  in  3-space.  Then  the  graph  of  the  upper  rim  h''{Q)  of  the  segments 
in  T  consists  of  at  most  0{n  a{n))  smooth  portions,  where  a(n)  is  the 
extremely  slowly  growing  inverse  Ackermann's  function.  Moreover,  h""  can 
be  calculated  in  0{n  a{n)  log  n)  time. 

Proof:  As  above  we  assume  that  no  segment  e,  is  such  that  its  horizontal  pro- 
jection e*  crosses  the  horizontal  ray  at  orientation  0  =  0  from  a*.  As  noted 
above,  this  involves  no  real  loss  of  generality  since  we  can  always  replace  t 
by  another  collection  having  at  most  2n  edges  which  satisfies  this  additional 
requirement  and  which  has  the  same  upper  rim  as  t.  Let  the  sequence  of 
smooth  portions  of  h'^  he  {hi,  .  .  .  ,h„),  arranged  in  counterclockwise  order 
along  S^  and  starting  at  0  =  0,  such  that  for  each  i^m  there  exists  a  segment 
e„_  €  T    such    that    all    the    rays    from    a    having    spherical    coordinates 
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(e,/2(e))  c  h^ 

pass  through  the  segment  e„, .  (If  in  such  a  smooth  portion  all  these  rays  pass 
through  two  segments  e,  e'  of  t,  which  can  happen  if  c,  e',  and  a  are  all  co- 
planar,  then  we  take  e„^  to  be  the  farther  of  these  two  segments.)  Let  U  be 
the  corresponding  sequence  (mi,«2.  •  •  •  .«m)-  We  claim  that  t/  is  an  (/j,3)- 
Davenport-Schinzel  sequence  in  the  terminology  of  [HS].  That  is,  U  satisfies 
the  following  properties 

(1)  Uii^Ui  +  i  for  each  i<m. 

(2)  For  each  pair  of  segment  indices  p  i^  q,  there  do  not  exist  five  indices 
ii<i2</3<i4</5  in  f/such  that  Ut^  =  Ut^=Ui^=p  and  Ui^=Ui^=q. 

Indeed,  the  first  property  is  obvious.  To  establish  the  second  property,  let 
Cpi^Cq  i  T  be  such  that  U  contains  an  alternation  of  the  form 
p  •  •  ■  q  •  '  •  p  •  •  •  q  •  •  •  p,  and  let  0i<92<e3  be  three  orientations  that 
correspond  to  the  three  middle  occurrences  q  •  •  •  p  •  •  •  q  in  this  subse- 
quence; in  particular,  the  rays  rej  =  (9i,/i(ei)),  rf,^  =  {B^,h{B2))  pass  through 
e^,  and  the  ray  re.  =  (02.''(02))  passes  through  ep.  It  is  easily  checked  that 
rej  and  r^^  both  pass  above  Cp  (in  the  strong  sense  that  the  projection  e* 
intersects  the  projections  rl^  and  rjj  of  these  rays,  and  that  above  each  of 
these  points  of  intersection  the  corresponding  ray  lies  above  Cp)  and  that  r^^ 
passes  above  e, .  But  then  it  follows  by  continuity  that  there  must  exist  two 
rays  r' ,  r"  (whose  projections  onto  the  j:->'  plane  have  orientations  lying  in 
(61,62)  and  in  (62,63)  respectively),  such  that  both  r'  and  r"  pass  through 
ep  and  through  e^.  This  implies  that  ep ,  c,  and  a  must  be  co-planar,  in  which 
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case,  by  our  convention,  only  one  of  Cp,  e^  should  be  used  to  represent  any 
of  the  smooth  portions  of  /z*"  in  which  it  passes  through  both  these  segments, 
a  contradiction  which  establishes  our  claim. 

U  is  thus  an  (n,  3)  Davenport-Schinzel  sequence,  and  by  the  results  of  [HS] 
the  length  of  U  is  at  most  0{na{n)).  The  techniques  in  [At],  [HS]  can  then 
be  used  to  calculate  h''  in  time  C>(na(n)  log  n),  using  a  straightforward 
divide-and-conquer  approach,  n 

Using  Proposition  1,  the  complexity  of  the  preprocessing  of  a  can  be 
estimated  as  follows.  We  note  that  if  t  is  a  collection  of  segments  which  is 
the  union  of  tj  and  t2  then  the  upper  rim  of  t  can  be  computed  by  "merging" 
the  upper  rims  of  tj  and  ti  in  time  linear  in  the  size  of  these  rims.  Thus  if 
we  calculate  the  upper  rim  of  the  entire  collection  t  of  the  edges  of  o- ,  so  that 
at  each  recursive  level  we  partition  the  current  set  of  edges  into  its  nearer 
half  and  further  half,  we  will  obtain,  in  total  time  O {n  a{n)  \og  n) ,  the  entire 
collection  of  upper  rims  that  we  need  to  store  in  our  search  tree  T.  (In  fact, 
only  half  of  the  rims  obtained  in  this  manner  need  be  stored  in  T\  see  Fig.  2.3 
for  an  illustration.)  The  total  size  of  all  these  upper  rims  is  easily  seen  to  be 
also  0{na{n)  log  n).  As  remarked  above,  computing  the  tree  T  that  sup- 
ports logarithmic-time  queries  requires  a  further  0{na.{n)  log  n)  time,  so  the 
total  time  used  to  build  T  is  0{n  a{n)  log  n).  In  summary,  we  have 

Theorem  2:  One  can  preprocess  a  polyhedral  terrain  cr  having  n  edges  in 
0{n  a(n)  log  n)  time  and  space,  so  as  to  support  C>(log  n)  ray  shooting 
queries  from  a  fixed  given  point. 
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Remark:  The  upper  bound  of  0{n  a(n))  on  the  complexity  of  the  upper  rim 
of  n  segments  is  tight  in  the  worst  case.  This  follows  from  a  recent  result  of 
Wiemik  and  Sharir  [WS]  that  the  upper  envelope  of  a  collection  of  n  seg- 
ments in  the  plane  can  consist  of  Cl{ncL{n))  subsegments.  Let  ei,  .  .  .  ,c„  be 
n  segments  whose  upper  envelope  consists  of  such  a  nonlinear  number  of 
subsegments,  and  which  are  all  drawn  in  the  upper  half  of  the  x—z  plane. 
We  can  construct  from  these  segments  a  polyhedral  terrain  as  follows.  Shift 
each  segment  e,  in  the  ^'-direction  to  the  plane  y  =  i,  and  associate  with  each 
shifted  segment  e*  a  sharp  wedge  having  e'  as  its  upper  edge.  The  terrain  tr 
is  then  defined  as  the  upper  envelope  of  all  these  wedges  and  of  the  x  —y 
plane;  see  Fig.  2.4.  Plainly  cr  has  0{n)  faces.  If  we  view  a  from  a  point  a 
lying  on  they-axis  sufficiently  far  away  from  these  segments,  then  it  is  easily 
checked  that  the  complexity  of  the  resulting  upper  rim  (of  the  entire  collec- 
tion of  segments  of  cr)  is  also  Vl{n  a(n)). 

3.  Visibility  of  a  polyhedral  surface  from  a  varying  viewing  point 

In  this  section  we  consider  an  extension  of  the  problem  studied  in  Section 
2,  in  which  the  point  a  from  which  the  surface  ct  is  viewed  is  allowed  to  vary. 
We  treat  here  only  the  simplest  extension  in  which  a  can  be  positioned  any- 
where along  a  fixed  vertical  line  L.  In  this  case  the  space  of  rays  that  need  be 
processed  is  3-dimensional,  and  is  conveniently  represented  by  LxS^,  so 
that  each  point  (fl,M)€Lx5^  represents  the  ray  emanating  from  a  in  direc- 
tion «  C  5^. 

The  fact  that  L  is  vertical  enables  us  to  obtain  a  preprocessing  algorithm 


-13- 

whose  time  complexity  is  0{n\4{n)\ogn);  it  produces  a  data  structure  of 
size  0{n\^(n))  which  supports  ray  shooting  queries  in  time  0(\og^n).  Our 
preprocessing  algorithm  actually  yields  an  explicit  representation  of  the  way 
in  which  a  can  be  viewed  from  L.  That  is,  it  produces  a  decomposition  of 
LXS^,  which  we  call  the  visibility  structure  of  ct  from  L,  into  maximal  con- 
nected cells,  so  that  with  each  cell  c  there  is  associated  a  face/^  of  ct  such  that 
all  rays  (a,u)  €  c  first  hit  cr  at  a  point  on/^.  We  show  that  in  the  worst  case 
the  boundaries  of  the  cells  in  such  a  decomposition  can  contain  fl(n^a(n)) 
vertices  and  edges,  so  that  our  algorithm  is  close  to  being  optimal  among 
algorithms  required  to  calculate  the  entire  visibility  structure.  There  remains 
of  course  an  intriguing  open  problem  whether  one  can  obtain,  as  in  Section  2, 
a  more  efficient  imphcit  representation  of  the  visibility  structure,  that  will 
still  support  fast  ray  shooting  queries.  We  also  show  that  if  L  is  not  vertical, 
there  are  cases  in  which  the  complexity  of  the  visibility  structure  of  LxS^ 
can  be  Cl{n^),  indicating  that  visibility  from  a  vertical  line  is  indeed  a  special 
favorable  case. 

Our  algorithm  proceeds  as  follows.  Since  all  points  a  6  L  have  the  same 
projection  a '  on  the  x  -y  plane,  the  order  among  the  edges  of  cr  as  defined 
in  Section  2  is  fixed  and  independent  of  the  choice  of  a  6  L.  Our  approach  is 
to  calculate  a  search  structure,  different  from  that  used  in  Section  2,  which  is 
defined  as  follows. 

Consider  the  collection  of  all  rays  in  Lx5^  whose  horizontal  projection 
has  some  fixed  orientation  Bq.  All  these  rays  lie  in  the  vertical  half  plane 
V(Qo)-  0  =  6o.    This  plane  intersects  ct  in  a  polygonal  line  rT(0o),  whose 
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c»mers  lie  on  edges  of  a  and  whose  edges  lie  on  faces  of  a.  Let  T  =  T(Qq) 
be  the  collection  of  (all  straight  segments  comprising)  the  shortest  paths  from 
a*  to  all  corners  of  cr(6o)  constrained  to  lie  in  V{Qq)  above  ct(0o)-  It  is  easily 
checked  (and  well  known,  as  a  special  case  of  shortest  path  trees  inside  sim- 
ple polygons  [GHLST])  that  the  edges  in  T  form  a  planar  tree  rooted  at  a*, 
that  there  are  0(n)  such  edges,  and  that  the  shortest  path  from  a*  to  each 
comer  v  of  cT(eo)  is  concave  (i.e.  all  its  bends  are  downwards).  We  represent 
points  in  V(eo)  by  the  remaining  pair  (r,z)  of  their  cylindrical  coordinates, 
with  7  being  the  vertical  axis.  See  Fig.  3.1  for  an  illustration  of  V  and  T. 

Consider  now  the  transformation  in  which  each  ray  emerging  from  a 
point  on  L  in  the  half -plane  V(eo)  is  mapped  into  a  dual  point,  as  defined  in 
Chazelle  and  Guibas  [CG2]  (namely,  each  ray  z  =  ar  +  b,  r>0,  is  mapped  to 
the  point  {—a,b)).  As  shown  in  [CG2],  in  this  dual  structure,  the  collection 
P{e)  of  all  such  rays  hitting  first  a  fixed  edge  e  of  (t(6o)  is  a  convex  polygon; 
its  vertices  are  the  duals  of  rays  passing  through  two  comers  of  ct(6o)  (one  of 
which  can  be  an  endpoint  of  e)  before  hitting  e.  Moreover,  the  special  struc- 
ture of  ct(0o)  implies  that  each  dual  polygon  P(e)  is  unbounded.  This  follows 
from  the  fact  that  if  a  ray  p  emerging  from  some  a  6  L  first  hits  cr(0o)  at  a 
point  q  ^  e,  then  all  rays  connecting  points  a'  i  L,  lying  above  a,  to  ^  have 
the  same  property.  Let  D(0o)  denote  the  partitioning  of  the  dual  plane  into 
these  polygons;  see  Fig.  3.2  for  an  illustration  of  this  stmcture.  Then  it  is 
easy  to  see  that  as  we  vary  0  slightly  near  0o,  the  combinatorial  structure  of 
the  partitioning  D{Q)  remains  unchanged,  unless  one  of  the  following  critical 
conditions  occur: 


-15- 

(i)  V(0o)  passes  through  a  vertex  of  a. 

(ii)  Two  adjacent  edges  in  T(Bq)  become  coUmear, 

Call  an  orientation  at  which  (i)  or  (ii)  occurs  a  critical  orientation. 
Qearly  there  are  only  0{n)  critical  orientations  of  the  first  type.  To  estimate 
the  number  of  critical  orientations  of  the  second  type  we  proceed  as  follows. 
For  each  edge  e  of  ct  and  for  each  orientation  6,  let  e{Q)  denote  the  intersec- 
tion of  e  with  V(0).  Fix  an  edge  c  of  ct  and  for  each  edge  ct  <  e  define  the 
function  /,(0)  to  be  the  azimuth  ^  of  the  ray  emerging  from  some  point  on  L 
and  passing  through  ci  and  through  the  point  e(6).  Define 
F,(e)  =  max  {  max //(e)  ,  <{>o(6)  },  where  <t>o(6)  is  the  azimuth  of  the  ray 


ei<e 


emerging  from  a'  towards  e{Q).  Note  that  if  F,{Q)  is  attained  by  some  // 
(i.e.  if  F,(e)  >  <})o(e)),  then  the  last  edge  in  the  shortest  path  from  a'  to 
e(9)  in  V(e)  (as  defined  above)  must  be  the  segment  e,(e)e(e).  Thus  6  is  a 
critical  orientation  of  the  second  type  at  which  the  last  edge  along  the  shor- 
test path  from  a*  to  e  (0)  changes  if  and  only  if  the  upper  envelope  F,  has  a 
"breakpoint"  at  6,  i.e.  F,  is  simultaneously  attained  by  two  functions/,,/^  (or 
//,  <t)o)  at  0.  For  each  edge  e  of  ct  let  1(e)  denote  the  angular  sector  spanned 
between  the  two  segments  connecting  a*  to  the  two  endpoints  of  c*.  Then  it 
is  easy  to  check  that  each  of  the  functions  /,(0)  is  partially  defined  and  con- 
tinuous over  the  connected  angular  interval  I{e)ni{ei)  (and  (t)o(0)  is 
defined  and  continuous  over  the  entire /(e)).  Moreover,  we  have 

Lemma  3:  Each  pair  of  functions //(0),  fj{Q)  intersect  in  at  most  two  points. 
Also  each  function /,(0)  intersects  4)o(0)  in  at  most  one  point. 
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Proof:  We  need  to  estimate  the  number  of  orientations  6  at  which  the  three 
points  c/(6),  ej{Q),  and  e(e)  are  collinear.  Consider  the  point  e(Q)  in  the 
standard  (r,r)  cylindrical  coordinate  frame  in  V(6).  Let  d  be  the  distance 
from  the  origin  to  the  line  \i  containing  the  projection  e*  ol  e.  It  is  straight- 
forward to  check  that  the  {r,z)  coordinates  of  c(e)  are 


(■ 


•y  cos  6  +  S  sin  9  X 


a  cos  9  -I-  p  sin  6     a  cos  9  -h  p  sin  9 
for  appropriate  constants  a,  P,  7,  8.    Similar  expressions  represent  e/(9), 

e^(9),  allowing  us  to  express  the  collinearity  of  these  three  points  by  an  equa- 
tion of  the  form 


DET 


d 7  cos  9  4-  S  sin  9 

a  cos  9  -I-  P  sin  9        a  cos  9  +  P  sin  9 
di  7/  cos  9-1-8/  sin  9 


a ;  cos  9  -I-  P;  sin  9      a^  cos  9  +  P/  sin  9 
dj  "ij  cos  9  4-  8y  sin  9 


=  0 


ttj  cos  9  4  Py  sin  9  ay  cos  9  4-  Py  sin  9 
It  is  easy  to  check  that  this  equation  can  be  rewritten  as  a  quadratic 
homogeneous  equation  in  cos  9,  sin  9,  or  rather  as  a  quadratic  equation  in 
tan  9.  We  thus  obtain  at  most  two  values  of  tan  9,  each  yielding  at  most  one 
value  of  9  in  the  angular  interval  1(e)  over  which  F^  is  defined  (whose  length 
is  clearly  less  than  tt).  This  implies  that  the  equation  //(9)  =  fj(Q)  has  at 
most  two  roots.  The  analysis  of  equations  of  the  form  /i(9)  =  <t)o(9)  is  the 
same,  except  that  the  last  row  of  the  above  determinant  is  [1,  0,  0]  (assuming 
that  a'  lies  at  the  origin).  This  leads  to  a  linear  homogeneous  equation  in 
cos  9,  sin  9  and  thus  yields  at  most  one  root  within  1(e).  □ 


As  shown  by  Atallah  [At],  the  total  number  of  breakpoints  along  the 
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upper  envelope  of  n  partial  functions  each  having  a  connected  graph,  such 
that  each  pair  of  them  intersect  in  at  most  s  points,  is  X,  +2(n).  Thus  the  max- 
imal number  of  breakpoints  in  F,  is  \^{n)  (which,  by  the  results  of  [Shi],  is 
0{n  a(n)*^(°  (")))).  Repeating  this  argimient  for  each  edge  e  of  a  we  con- 
clude that  the  total  nimiber  of  critical  orientations  of  the  second  type  is 
0{n\4{n)).  Furthermore,  the  divide-and-conquer  technique  in  [At]  can  be 
used  to  obtain  all  these  critical  orientations  in  overall  time  0{nX^{n)  log  n). 
Actually,  below  we  will  use  an  alternative  technique  to  calculate  these  orien- 
tations within  the  same  time  bound. 

We  turn  to  describe  the  data  structure  that  will  be  used  to  represent  the 
visibility  structure  we  seek.  First,  we  define  a  total  order  on  the  faces 
S\Ji,  •  •  ■  ,/n  of  a  so  that  if  i  <;,  then  any  ray  enemating  from  L  and  inter- 
secting both//  and^,  intersects//  first.  Such  an  order  exists  if  we  make  the 
assumption,  similar  to  that  made  in  the  previous  section,  that  no  face  of  cr  is 
cut  by  the  vertical  half-plane  6  =  0.  (To  enforce  this  property,  we  simply 
split  each  face  that  is  cut  by  that  half-plane  into  two  subfaces.)  We  can  then 
compute  this  order  in  linear  time  as  follows.  Project  the  faces  of  ct  onto  the 
x—y  plane.  Define  a  partial  order  between  the  faces  so  that  fi<fj\f  they  are 
adjacent  and  all  horizontal  rays  from  a  *  that  meet  the  edge  between  /J  and  fj 
intersect  ft  before  they  intersect  fj.  Then  complete  this  order  to  a  total 
order  by  topological  sorting  in  0{n)  time. 

Next,  we  need  to  study  the  dual  structure  D  (6)  more  carefully.  We  will 
use  {r,2)  coordinates  for  the  primal  half-plane  V(e),  and  {r,z)  coordinates 
for  the  dual  structure  Z)  (0).   Consider  the  faces  /,  of  ct  which  are  crossed  by 
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the  half -plane  V(0),  and  let  e/  denote  the  intersection  of  //  with  that  half- 
plane.  To  each  such  segment  et  there  corresponds  a  polygon,  which  we 
denote  by  Pi,  in  D  (6).  As  noted  above,  Pi  is  an  unbounded  convex  polygon, 
and  the  unbounded  rays  that  it  contains  all  have  positive  slope.  It  is  also 
easily  checked  that  the  order  of  the  polygons  in  D  (6),  as  they  intersect  a  suf- 
ficiently distant  horizontal  line  from  right  to  left,  is  the  same  as  the  order  of 
the  corresponding  faces  /,  of  ct.  In  addition,  every  bounding  edge  (of  a 
polygon  in  D  (9))  consists  of  dual  points  of  rays  that  pass  through  a  vertex  of 
(T(e),  so  that,  arguing  as  above,  each  such  edge  must  have  positive  slope. 
Thus  the  line  L  that  separates  regions  P\,P2,  •  •  •  >^m  from  regions 
P„  +  i,  .  .  .  ,P„  is  monotone  with  respect  to  both  the  z  and  the  f-axes. 

We  store  D(Q)  in  a  data  structure  that  supports  planar  point  location 
queries  in  C>(log^n)  time,  as  follows.  (There  are  of  course  many  well-known 
optimal  data  structures  that  will  support  O(log^z)  time  point  location  queries, 
but  it  does  not  appear  possible  to  modify  these  data  structures  to  allow  for 
varying  6  without  sacrificing  their  space  efficiency.)  We  create  a  complete 
binary  tree  Q  with  n  leaves,  one  for  each  of  the  regions  Pi,  .  .  .  ,P„.  We  call 
Q  the  region  tree.  The  symmetric  order  of  the  regions  in  Q  is  the  same  as  the 
total  order  of  the  corresponding  faces  of  a.  (Since  not  all  faces  of  cr  neces- 
sarily intersect  V(0),  some  of  the  regions  P,  can  be  empty  for  a  particular  6.) 
At  each  internal  node  v  of  G  we  store  the  polygonal  line  separating  the 
regions  stored  in  the  left  subtree  of  v  from  the  regions  stored  in  its  right  sub- 
tree. More  precisely,  this  line  consists  of  those  segments  which  separate  a 
region  stored  in  the  left  subtree  of  v  from  a  region  stored  in  the  right 
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subtree.  Consequently,  this  separating  line  may  consist  of  several  uncon- 
nected segments.  (For  example,  consider  the  line  separating  Pi,  .  .  .  ,P„/4 
from  P„/4+i,  .  .  .  ,Pni2-  It  might  be  that  at  certain  points  the  regions 
Pi,  .  .  .  ,P„/4  are  adjacent  to  the  line  separating  Pi,  .  .  .  ,P„/2  from 
P„/2+i,  •  •  •  ,Pn-  For  such  points  there  is  no  line  separating  Pi,  .  .  .  ,P„/4 
from  P„/4+i,  .  .  .  ,P„/2-)  Thus  each  edge  in  D(e)  is  stored  in  exactly  one 
separating  line;  specifically,  suppose  edge  e  separates  P,  and  Pj\  let  i:^k-2'<j 
be  the  integer  of  this  form  with  the  largest  possible  r;  then  e  belongs  to  the 
line  separating  P(t_  1)2'+ 1,  .  .  .  ,Pkr  fromPt.2'+i,  .  .  .  ,P(jk  +  i)2^-  SinceD(e) 
contains  0{n)  edges,  the  overall  complexity  of  all  these  separating  lines  is 
also  0{n).  Fig.  3.3  shows  the  collection  of  separating  lines  for  the  structure 
D(0)  shown  in  Fig.  3.2.  Each  separating  line  is  stored  in  a  balanced  binary 
search  tree,  ordered  by  f  (or  2)-coordinate. 

We  now  explain  how  to  perform  point  location  in  D(e).  Suppose  that 
we  are  locating  point  {a,b).  We  start  by  determining  which  side  of  L,  the 
separating  line  at  the  root  of  Q,  contains  {a,b).  We  do  this  by  determining 
the  edge  e  in  L  that  contains  a  point  with  f-coordinate  a,  by  means  of  a 
binary  search.  We  then  determine  for  {a,b)  on  which  side  of  e  it  lies,  and 
recursively  search  the  left  subtree  if  {a,  b)  is  below  L,  or  the  right  subtree  if 
it  is  above  L.  In  general,  there  may  be  no  such  edge  e\  but  then  we  know 
that  (a,  b)  is  on  the  same  side  of  L  and  of  J,  where  J  is  the  separating  line 
stored  at  the  parent  of  the  current  node.  Thus  we  can  perform  the  planar 
point  location  in  D  (9)  in  overall  O(log^n)  time. 
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However,  we  need  to  produce  a  data  structure  which  supports  fast  point 
location  queries  in  D  (6)  for  arbitrary  0.  Actually,  the  data  structure  we  have 
just  described  is  appropriate  for  all  6  in  the  range  (Gi,  62)  where  61  and  62 
are  two  adjacent  critical  orientations,  as  defined  earlier.  Moreover,  when  we 
cross  a  critical  orientation  6  of  the  second  type,  the  structure  Z)  (6)  undergoes 
only  0(1)  changes,  consisting  of  addition  and  deletion  of  edges;  if  6  is  one  of 
the  0(n)  type  (i)  critical  orientations,  the  total  change  in  D{B)  is  at  most 
0(n).  This  easily  implies  that  the  total  number  of  changes  in  D{Q),  and  in 
the  collection  of  all  separating  lines,  is  at  most  0{n\^{n)). 

At  each  node  v  of  the  actual  region  tree  Q,  we  store  the  entire  collection 
of  separating  lines  for  v,  as  they  vary  as  6  increases  from  0  to  2tt,  using  the 
similar  lists  data  structure  of  [Co]  or  [DSST].  The  space  requirement  for  the 
entire  tree  Q  is  thus  0(n  +  71X4(11)),  and  Q  now  supports  binary  search  over 
each  separating  line  in  time  0(log/i),  yielding  overall  0{\og^n)  time  for  a 
point  location  query,  for  arbitrary  6. 

It  remains  to  describe  how  to  build  the  modified  data  structure  Q.  This 
is  readily  accomplished  by  sweeping  the  half-plane  V(9)  from  6  =  0  to  It:, 
maintaining  the  tree  7(6)  and  the  dual  structure  D(e),  and  updating  the  tree 
Q  at  each  value  of  6  in  which  7(6)  and  D  (6)  change.  The  initial  value  of  Q 
is  the  search  structure  for  the  map  D  (0) . 

The  shortest  path  tree  7(6)  is  stored  in  a  manner  consistent  with  its 
definition;  that  is,  the  parent  of  each  comer  v  of  or(6)  is  its  predecessor  along 
the  shortest  path  from  a*  to  v.  The  preorder  of  the  nodes  of  7(9)  coincides 
with  their  r-order  as  corners  of  V(e).  As  6  increases,  7(6)  can  change  in  two 
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ways:  (i)  Several  nodes  are  removed  and  replaced  by  other  nodes  (this  hap- 
pens at  critical  orientations  of  type  (i)).  (ii)  The  parent  of  a  node  v  is 
changed  to  the  nearest  sibling  on  the  left  of  v  or  to  its  grandparent,  with  the 
second  type  of  change  occurring  only  if  v  is  the  rightmost  child  of  its  parent 
(this  happens  at  critical  orientations  of  type  (ii)).  Since  there  are  only  0{n) 
changes  of  type  (i),  we  simply  re-calculate  T(Q),  in  0(n)  time,  at  each  of  the 
type  (i)  critical  orientations.  Each  of  the  updates  of  type  (ii)  is  easily  accom- 
plished in  0(1)  time.  Thus  the  total  time  for  all  updates  of  7(6)  is 
0(n\,in)). 

To  perform  these  updates  during  our  sweep,  we  need  to  know  which  is 
the  next  update  to  occur.  For  this,  we  maintain  a  heap  of  candidate  critical 
orientations.  In  this  heap  we  store  the  0{n)  orientations  defined  by  the  0{n) 
triples  of  corners  of  cr(6)  that  might  cause  a  change  of  type  (ii),  namely:  for 
each  pair  of  adjacent  siblings  in  r(6)  we  store  the  orientation  6  at  which  the 
rightmost  of  the  two  siblings  becomes  a  child  of  its  sibling,  rather  than  of  its 
parent  (i.e  the  orientation  at  which  the  three  comers  defining  the  parent  and 
the  two  siblings  become  colinear);  similarly  we  store  the  orientation  at  which 
a  triple  comprising  a  rightmost  sibling,  its  parent  and  grandparent  become 
colinear.  (For  nodes  that  are  children  of  the  root  of  r(0)  the  second  type  of 
change  cannot  occur,  so  we  do  not  have  to  store  any  corresponding  orienta- 
tion.) We  also  add  to  the  heap  the  0{n)  critical  orientations  at  which  changes 
of  type  (i)  occur.  Note  that  at  each  critical  orientation  of  type  (ii)  there  are 
0(1)  updates  (insertions  and  deletions)  to  be  made  to  the  heap;  they  require 
C>(log  n)  time.   At  each  orientation  of  type  (i)  we  can  simply  re-calculate  the 


-22- 

heap  in  time  0{n  log  n).  So  overall,  maintaining  the  heap  takes  time 
0{n\4{n)  log  n). 

Finally,  we  explain  how  to  construct  the  point  location  data  structure  Q. 
As  observed  above,  each  change  to  7(6)  produces  corresponding  changes  to 
D(0).  Each  of  these  changes  consists  of  the  addition  or  removal  of  a  single 
known  edge  from  one  of  the  separating  lines.  Furthermore,  if  e  is  an  edge 
separating  two  regions  P,  and  Pj ,  which  is  being  added  to  (or  removed  from) 
D(e),  we  can  calculate  from  the  indices  /  and  j,  in  the  manner  explained 
above,  which  node  v  of  Q  should  store  (or  already  stores)  e.  It  remains  to 
insert  (or  delete)  e  into  (from)  the  separating  line  L  stored  at  v.  But  this  is  a 
standard  operation  for  the  similar  lists  data  structure  [Co],  [DSST],  and  can 
thus  be  done  in  0(log  n)  time.  Thus  the  data  structure  Q  can  be  built  in 
0{n  X4(n)  log  n)  time. 

To  sum  up,  we  have  shown 

Theorem  4:  The  visibility  structure  of  a  polyhedral  terrain  ct  with  n  faces  in 
Lx  S^,  for  a  given  vertical  line  L,  is  of  size  at  most  0(nX4(n)).  It  can  be  cal- 
culated in  time  0{nX^{n)  log  n)  and  stored  into  a  data  structure  of  size 
0{n'K4{n))  which  supports  ray  shooting  queries  in  C>(log^/i)  time. 

Example  5:  We  do  not  know  whether  the  bound  0{n\4(n))  is  tight  in  the 
worst  case,  but  the  following  example  shows  that  the  complexity  of  the  visi- 
bility structure  can  be  n(n'K^{n))  =  Ci{n^a{n)).  In  this  example  we  first 
take,  as  in  Section  2  (and  following  [WS]),  k  segments  whose  upper  rim  as 
viewed  from  a*  consists  of  Cl{ka{k))  subsegments.    Each  of  these  segments 
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lie  in  a  vertical  plane,  and  these  planes  are  all  parallel  and  sufficiently  close  to 
one  another  so  that  the  combinatorial  structure  of  their  upper  rim  remains 
unchanged  as  a  varies  along  L  from  a*  to  a  sufficiently  high  point  a". 
Moreover,  by  scaling  down  the  y-coordinate  in  the  construction  of  [WS],  we 
can  assume  that  all  the  breakpoints  in  the  upper  rim  of  these  segments 
correspond  to  rays  passing  sufficiently  near  some  fixed  horizontal  line  M. 
We  can  now  complete  ct  by  it  additional  long  horizontal  parallel  segments 
lying  beyond  the  first  it  edges  such  that  the  faces  they  bound  form  a  convex 
surface  (bending  upwards),  and  such  that  as  a  varies  along  L  from  a*  to  a**, 
the  rays  shooting  from  a  and  passing  through  M  will  encounter  each  of  the 
latter  it  edges.  See  Fig.  3.4  for  an  illustration  of  such  a  terrain.  It  is  now 
clear  that  for  such  a  terrain  ct,  the  complexity  of  the  visibility  structure  is 
Cl{n'^a.{n)).  In  particular,  it  contains  n(n^a(n))  vertices,  i.e.  rays  that  pass 
through  three  edges  of  ct  simultaneously. 

Remark:  The  requirement  that  the  viewing  line  L  be  vertical  can  be  trivially 
relaxed  by  requiring  only  that  each  line  parallel  to  L  intersects  ct  in  precisely 
one  point.  In  other  words,  the  above  analysis  will  continue  to  apply  for  non- 
vertical  lines  L,  provided  that  when  the  coordinate  system  is  rotated  so  as  to 
make  L  vertical,  ct  remains  the  graph  of  a  single  valued  polyhedral  function. 

Example  6:  We  conclude  this  section  by  an  example  in  which  the  line  L  is  not 
vertical  and  the  complexity  of  the  visibility  structure  in  Lx5^  is  ft(«^).  Let  L 
be  a  horizontal  line,  and  let  ct  consist  of  two  parallel  rows  of  it  sharp  peaks 
each,  beyond  which  k  long  horizontal  parallel  edges  are  placed  so  that  the 
faces  they  bound  form  a  convex  surface  (bending  upwards).  L  is  placed  in 
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front  of  the  first  row  of  peaks  and  parallel  to  it.  See  Fig.  3.5  for  an  illustra- 
tion of  this  construction.  Now  for  each  pair  of  peaks  p,  p' ,  one  lying  in  the 
first  row  of  peaks  and  the  other  in  the  second  row,  there  exists  a  point  a  i  L 
which  is  co-planar  with  two  respective  parallel  sides  e,  e'  of  these  peaks. 
Qearly  in  this  plane  there  exist  it  rays  shooting  from  a,  each  of  which  passes 
through  e,  e'  and  one  of  the  further  it  horizontal  edges  of  cr.  Thus  the  total 
number  of  rays  shooting  from  L  and  passing  through  three  edges  of  cr  simul- 
taneously is  n(it-')  =  n(/i^),  as  claimed.  (We  are  grateful  to  Ady  Wiemik 
for  pointing  out  this  example.)  Of  course,  this  does  not  preclude  the  possibil- 
ity of  an  implicit  and  more  efficient  representation  of  LX-S^  which  still  sup- 
ports fast  ray  shooting  queries. 

4.  Visibility  of  an  entire  terrain  from  several  points. 

We  conclude  the  paper  with  a  few  remarks  on  the  problem  of  posting 
several  viewing  stations  on  the  surface  of  a  polyhedral  terrain  a  with  n  faces 
so  that  each  point  on  cr  is  visible  from  at  least  one  of  these  stations.  We 
show  that  for  the  case  of  a  single  station,  its  existence  (and  location)  can  be 
determined  in  time  0{n  log  n)  in  a  straightforward  manner.  On  the  other 
hand,  the  problem  of  determining  the  smallest  number  of  stations  that  cover 
the  entire  terrain  is  A^P-hard. 

Consider  first  the  case  of  a  single  viewing  station  a  €  ct.  For  each  face  /, 
of  cr  let  Ht  be  the  (closed)  half-space  bounded  by  the  plane  containing  //  and 
lying  above  it.  We  claim  that  ct  is  entirely  visible  from  a  if  and  only  if 

n 

a  ^   n  Hi.  Indeed,  for  each  face//,  if  a  can  see  points  on  //  it  must  clearly  lie 
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Hi.  Conversely,  let  a  ^   fl  ///  and  suppose  that  some  point  z  on  a  face  fi  of  cr 

is  not  visible  from  a.  Then  the  ray  az  must  intersect  cr.  Let  b  be  the  last 
intersection  of  az  with  cr  before  reaching  z,  and  suppose  b  i  fj.  Then  it  is 
easily  seen  that  a  is  not  in  //y,  a  contradiction  that  proves  our  claim. 

It  thus  follows  that  there  exists  a  viewing  station  that  sees  the  entire  sur- 
face cr  if  and  only  if  6-  =  cr  fl  (fl  Hi)  is  nonempty,  and  then  any  point  in  cr 
can  serve  as  such  a  viewing  station. 

To  check  whether  a  is  non-empty,  and  also  to  find  a  point  on  it,  we  can 

n 

proceed  as  follows.  First  calculate  K  =   fl  ///  in  time  0{n  log  n)  using  the 

algorithm  in  [PM].  Suppose  that  cr  is  non-empty  and  that  it  contains  a  point  a 
lying  in  the  interior  of  some  face  /,  of  cr.  Let  g,  be  the  face  of  K  co-planar 
with  //  {gt  may  degenerate  to  an  edge  or  a  vertex  of  K).  Then  gt  must  be 
contained  in  //.  Indeed,  if  gi  contained  a  point  b  (.  fi  then  the  segment  ab 
must  lie  entirely  on  or  above  cr  and  must  also  cross  an  edge  (or  vertex)  of  // . 
But  then  the  face  of  a  adjacent  to  //  across  that  edge  (or  some  face  of  cr  adja- 
cent to  that  vertex)  will  not  be  visible  from  a.  Similarly,  if  &  contains  a  point 
a  lying  on  an  edge  e,  of  ct  then  the  edge  of  K  collinear  with  ei  must  be  con- 
tained in  €{.  This  argument  proves  that  d  is  non-empty  if  and  only  if  cr  con- 
tains a  vertex  of  K.  This  latter  property  is  easy  to  check  using  a  standard 
sweeping  technique  applied  to  the  x  —y  projection  ct*  of  a  and  to  the 
corresponding    projections    of    the    vertices    of    K,    which    runs    in    time 
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0{n  log  n).  Thus  we  have 

Proposition  7:  Given  a  polyhedral  terrain  ct  with  n  faces,  one  can,  in  time 
0{n  log  n),  detennine  whether  ct  contains  a  point  which  can  see  the  entire 
surface  ct,  and  find  such  a  point  when  it  exists. 

Consider  next  the  problem  of  determining  the  smallest  number  of  view- 
ing stations  on  a  needed  to  cover  the  entire  surface  of  o-  (in  the  sense  defined 
above).  We  show  that  this  problem  is  NP-haid  by  reduction  from  Satisfiabil- 
ity. 

Let  F  be  a  CNF  formula  with  clauses  Ci,  .  .  .  ,C„,  and  variables 
xi,  .  .  .  ,x„.  We  reduce  the  satisfiability  problem  for  F  to  the  problem  of 
determining  whether  a  certain  terrain  cr  with  0{mn)  faces  can  be  completely 
viewed  by  mn  points  on  it. 

The  surface  ct  to  be  viewed  is  a  horizontal  plane  on  which  certain  walls 
are  erected  and  certain  pits  are  dug  out.  There  are  n  — 1  parallel  walls;  they 
separate  n  rows,  one  per  variable.  The  pits  occur  in  the  rows;  There  are  2m 
pits  per  row  arranged  in  a  circular  fashion  as  shown  in  Fig.  4.1,  which 
displays  the  view  of  a  row  from  above.  The  upper  rims  of  the  pits  are  all 
quadrilateral,  and  the  rims  of  each  pair  of  adjacent  pits  in  the  same  row  have 
a  common  vertex,  which  we  call  a  peak.  We  create  m  columns,  perpendicular 
to  the  rows,  one  per  clause.  Each  column  will  cut  through  all  the  walls.  The 
cross-section  of  a  pit  lengthwise  (along  the  row)  is  shown  in  Fig.  4.2. 

The  idea  is  that  only  from  a  peak  are  the  interiors  of  the  two  incident 
pits  completely  visible,  and  only  from  the  boundary  or  interior  of  a  pit  can 
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we  see  the  whole  pit.  This  is  achieved  by  making  the  pit  very  deep  relative 
to  the  height  of  the  walls  (e.g.  as  shown  in  Fig.  4.3).  Thus  to  see  the  entire 
bottom  of  all  the  2m  pits  in  a  row,  at  least  m  viewing  points  will  always  be 
needed.  Moreover,  let  pi,  .  .  .  ,p2m  be  the  peaks  in  some  row  r.  Then  to  be 
able  to  view  all  the  pits  of  r  with  exactly  m  points,  they  will  have  to  be  placed 
at  every  other  peak  (i.e.  either  at  the  odd-numbered  peaks  or  at  the  even- 
numbered  peaks  of  that  row).  Assuming  row  r  corresponds  to  variable  x^, 
the  choice  of  even  peaks  for  the  viewing  points  in  r  will  correspond  to  setting 
Xr  =  true,  while  the  choice  of  odd  peaks  for  these  points  will  correspond  to 
setting  X  ==  false.  If  clause  Cj  of  our  formula  does  not  contain  x^  or  Xr  then 
the  peaks  P2j-i,  Pij  in  row  r  are  placed  outside  of  any  column;  if  Cj  contains 
Xr  (resp.  Xr)  then  peak  p2j  (resp.  P2j-i)  is  placed  inside  the  y-th  column 
while  P2J-1  (resp.  p2j)  is  placed  outside  of  any  column.  All  this  can  be 
easily  arranged  by  appropriately  varying  the  lengths  of  the  pits  in  the  r-th 
row. 

A  wall  has  a  triangular  cross-section  (parallel  to  the  columns).  Thus  a 
side  of  a  wall  can  be  viewed  in  its  entirety  from  any  peak  in  either  of  the 
adjacent  rows.  Thus  no  additional  visibility  points  are  required  for  viewing 
the  sides  of  the  walls  parallel  to  the  rows.  The  cross-section  of  a  wall  (paral- 
lel to  a  column)  is  shown  in  Fig.  4.4. 

Each  column  makes  a  vertical  cut  through  each  of  the  walls.  The  cross 
section  along  one  wall  (parallel  to  a  row)  showing  the  columns  is  given  in 
Fig.  4.5.  Finally,  we  ensure  that  each  column  is  completely  visible  from  any 
peak  occurring  in  the  column,  but  not  from  any  collection  of  peaks  outside 
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the  column.  This  is  achieved  by  making  the  walls  wide  and  the  columns  nar- 
row (the  first  and  last  columns  do  not  need  to  be  narrow,  and  in  fact  are 
bounded  on  just  one  side).  A  view  from  above  of  the  resulting  terrain  for 
the  formula 

F  =  (xi  Va:2  VJ3)  A(Ji  Vj:3  VJ4) 
is  shown  in  Fig.  4.6. 

Thus  the  reduction  results  in  a  polyhedral  terrain  cr  having  0{mn)  faces 
which  requires  at  least  mn  viewing  points  to  see  all  of  it.  Moreover,  if 
exactly  this  number  of  viewing  points  is  used,  all  of  them  must  be  placed  at 
peaks  (in  order  to  see  all  the  pits)  in  the  manner  described  above.  Further 
we  must  have  at  least  one  viewing  point  per  column.  But  this  can  occur  (with 
only  mn  viewing  points)  if  and  only  if  the  formula  F  is  satisfiable. 

We  have  thus  shown  that  Satisfiability  is  reducible  to  the  problem  of 
determining  whether  a  polyhedral  terrain  a  can  be  viewed  in  its  entirety  from 
a  given  number  of  viewing  points;  we  can  therefore  conclude 

Theorem  8:  It  is  A^P-hard  to  determine,  for  a  given  polyhedral  terrain  cr  with 
n  faces  and  a  given  integer  k,  whether  there  exist  k  viewing  stations  on  a 
which  collectively  see  the  entire  surface  of  cr. 

Remark:  Our  reduction  actually  shows  that  it  is  A^P-complete  to  determine, 
for  a  polyhderal  terrain  cr  and  an  integer  k,  whether  there  exist  k  vertices  of 
a  which  collectively  see  the  entire  surface  cr.  This  latter  problem  is  in  NP 
because,  after  having  guessed  the  it  vertices,  we  can  compute  in  polynomial 
time  the  portion  of  ct  visible  from  each  of  these  vertices,  and  then  check  that 
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the  union  of  these  portions  is  the  entire  surface  cr. 

5.   Conclusion 

In  this  paper  we  have  analyzed  a  collection  of  3-D  visibility  problems 
involving  viewing  a  polehedral  terrain  from  a  point  or  from  a  collection  of 
points.  As  far  as  we  know,  not  much  work  has  been  done  so  far  concerning 
detailed  analysis  of  the  computational  complexity  of  these  problems,  and  we 
hope  that  the  results  presented  here  will  spur  further  progress  in  this  area, 
which  we  believe  to  have  many  applications  in  computer  graphics,  geographi- 
cal data  processing,  and  other  related  image  processing  problems. 

Among  the  immediate  open  problems  that  are  suggested  by  our  results 
are: 

(1)  Preprocess  an  arbitrary  polyhedral  surface  (not  necessarily  a  terrain),  or 
even  an  arbitrary  collection  of  disjoint  polyhedral  surfaces,  so  as  to  support 
fast  ray-shooting  queries  from  a  fixed  viewing  point.  Can  this  be  done  with 
efficiency  comparable  to  that  obtained  in  Section  2? 

(2)  Obtain  sharp  bounds  for  the  combinatorial  complexity  of  the  visibility 
structure  for  a  terrain  cr  with  n  faces,  when  the  viewing  point  a  can  vary  in 
more  general  subspaces  L  of  E^.  In  particular  what  is  the  complexity  of  this 
structure  when  L  is  the  entire  half -space  lying  above  ct,  or  when  L  is  a  verti- 
cal half-plane  lying  above  ct? 

(3)  Can  the  results  of  Section  3  be  improved  to  obtain  implicit  and  more 
space-efficient  representation  of  the  visibility  structure  in  L  x  5^  for  a  verti- 
cal line  L? 
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(4)  Recently,  Plantinga  and  I>yer  [PD]  have  considered  similar  visibility  prob- 
lems concerning  parallel  views  of  an  arbitrary  polyhedral  scene  from  points  at 
infinity,  giving  a  worst-case  0(n'^)  algorithm  for  the  calculation  of  the  result- 
ing visibility  structure.  This  suggests  several  open  problems,  including  the 
problem  of  obtaining  sharp  bounds  on  the  complexity  of  the  visibility  struc- 
ture for  a  polyhedral  terrain  from  points  at  infinity,  and  the  problem  of  find- 
ing implicit  and  more  space  efficient  representations  of  this  structure  (for  ter- 
rains as  well  as  for  arbitrary  polyhedral  scenes). 
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Figure  1.1.  Visibility  of  a  terrain  from  a  point  with  fi(n^)  complexity. 
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(b) 


Figure  2.1(a).  The  order  of  edges  Figure  2.1(b).  A  possible  circularity. 


Figure  2.2.  The  upper  rim  of  a  terrain. 
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Figure  2.3.  The  tree  of  upper  rims  for  a  terrain  with  16  edges. 


Figure  2.4.  The  lower  bound  construction  for  upper  rims. 


y 


V(0n) 


y 


X 


/ 


/^     /" 


V(ao) 


^ 


^ 


\ 


\ 


©=er 


Figure  3.1.  V(eo)  and  r(eo). 
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Figure  3.2.  The  structure  D  (Bq)  for  the  structure  of  Figure  3.1. 
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Itie  internal  nodes  of  Q,  each  labeled  by  i:i+l,  where  P^  is 
the  largest-numbered  region  in  the  left  subtree  of  the  node. 


Node  of  Q 


Separating  edges  stored  at  the  node 


(i,j)  designates  the  edge  separating  P^  and  P. 
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Figure  3.3.  The  region  tree  Q  and  the  separating  lines. 


Figure  3.4.  The  lower  bound  construction  for  visibility  from  a  vertical  line. 
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Figure  3.5.  The  lower  bound  construction  for  visibility  from  a  horizontal  line. 
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Figure  4.1.  Schematic  arrangement  of  pits  in  a  row. 
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Figure  4.2.  A  vertical  cross-section  of  a  pit. 
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Figure  4.3.  A  vertical  cross-section  across  a  row. 


Figure  4.4.  A  vertical  cross-section  of  a  wall. 
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Figure  4.5.  A  wall  with  its  column  appertures  viewed  from  an  adjacent  row. 


Figure  4.6.  A  view  from  above  of  the  terrain  constructed  for 

F  =  (jTi  V  ^2  V  J3)  A  (Ji  V  j:3  V  J4). 
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